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Abstract 


A coordinate-free definition of complex-structure multidimensional item response theory 
(MIRT) for dichotomously scored items is presented. The point of view taken emphasizes the 
possibilities and subtleties of understanding MIRT as a multidimensional extension of the classical 
unidimensional item response theory models. The main theorem of the paper is that every 
monotonic MIRT model looks the same; they are all trivial extensions of univariate item response 
theory. 
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1 Introduction 


Complex-structure multidimensional item response theory (MIRT) is built on the idea that 
a single item, however simple it might be, carries the possibility of an inner structure. In usual 
terminology one speculates that it is possible to measure several cognitive areas with one item. 
The number of cognitive areas so measured may vary among items, even though usual models 
assume that it is fixed for a collection of items (a test) and let a factor analysis type procedure 
decide on the number and mixture of cognitive areas measurable by the items. 

The point of view taken in this note is that any uni dimensional item response theory (IRT) 
model can be thought of as a specialization of an MIRT model. Hence, the major task is to identify 
how much of the well-established tools and nomenclature of unidimensional IRT can be preserved 
in the multidimensional context and, from the other direction, how different multidimensional 
notions may specialize to the same unidimensional entity. When the latter happens, that is 
when two different multidimensional objects yield the same unidimensional specialization, then 
both multidimensional notions could be considered proper generalizations of the underlying 
unidimensional quantity. A careful study should then be devised to decide which generalization is 
more appropriate with respect to the application at hand. 

There is, on the other hand, the possibility of not finding proper multidimensional 
generalization for some unidimensional notions. This topic also deserves careful research and 
understanding. 

This paper only considers what is termed complex-structure MIRT. Usually, IRT models have 
two components: the item likelihood and the population distribution. In simple-structure MIRT, 
one item represents only one dimension and without a multivariate population distribution the 
entire likelihood of the model would factor as a product of univariate pieces. In complex-structure 
MIRT, this factorization is impossible, by definition, irrespective of population model chosen. 

The structure of the paper is as follows: a short overview of unidimensional IRT is followed 
by the absolute, that is, coordinate-free definition of MIRT. The connection with the usual 
approach is also shown via a discussion of two widely accepted models. Then the development 
of the main thesis follows. This paper proves that MIRT models are all alike and they all can 
be obtained as a trivial extension of an appropriate unidimensional IRT model. Two sections on 
some thoughts about capturing cognitive dimensions and on understanding the role of the notion 
of dimension-wise independence close the presentation. 
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2 Unidimensional Item Response Theory 

To make the generalization to the multidimensional framework easier, some features of 
unidimensional IRT are summarized first. Measurement takes place during the formation of 
the response matrix X £ Af/v x /(N) with elements x' n * £ N for student n = 1 ,,N and item 
i = 1,..., I. In a dichotomous setting (which is assumed throughout the paper to simplify the 
presentation), x n i = 1 if student n responded correctly to item i , otherwise it is zero. As a major 
simplification of the modeling of the cognitive process, it is assumed that students response to an 
item is stochastically determined by students ability 9 n and item parameters Pi := (aj,6j,Cj) via 
the item response function (Birnbaum, 1968). 

P 3p \0 n , Pi) := Prob (x ni = 1 | Q n , Pi) = c* + 1 % , . ■ (1) 

—L g Uiyun ) 

There are, of course, many different item response functions in use: the three-parameter logistic 
(3PL) model is chosen here only as an illustration. The other substantial simplification used in 
building the model is the assumption of independence of conditional probabilities P^f 1 across an 
arbitrary subset S C {1,..., N} x {1,..., 1} of student-item pairs. 

The two most popular models built out of these blocks are the joint unidimensional IRT and 
the marginal unidimensional IRT. Joint IRT states that the total likelihood depends explicitly on 
the ability of the given students: 

Li° M (X-e,B) = []i 53pl («n,ft) I ”(l - P Zv \e n ,Pi)Y~ Xni (2) 

n,i 

with corresponding log-likelihood: 

£j° int (X;0,6) = Y,Xnilog(P 3p \e n ,Pi)) + (1 - Xni)log(l - P 3p \0 n ,Pi)). (3) 

n,i 

Here, 0 and B are the collections of all abilities and item parameters, respectively. 

In marginal theory, the likelihood depends only on the distributional properties of student’s 
population: 

L marg( X ;£ i$) = TT f TT p3pl(^ £)*„, (1 _ p3pl(^ (0) (4) 

n J R i 

with log-likelihood 

£ mar s(X;H,<k) = Viog / TT P 3p \e,pp x -pi-p 3p \e,Pi)) 1 ~ x -dii n (d) 1 (5) 

i 
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where n n is the density measure of student n over M and $ is the collection of distributional 
parameters for student’s ability. In a parametric setting, usually fi n is given as 
dnn{0) = <p n (0) d# with some distribution function ip n . 

The quantities 

L*(X-,6,B) = Y[P 3p \0, Pi) Xni {l - P 3p \0, fa)) 1 -*™ (6) 

i 

and 

Lf(X;Q,P) = l[P 3 v\e n ,py^(l - P 3 v\e n ,P)y- x ^ (7) 

n 

are the student and item likelihoods, respectively. 

It is worthwhile to analyze the shape of the student likelihood function. It is a product of the 
conditional probabilities of the actual responses over all the items administered to the student. 
As a function of 0, the probability of the correct response increases when the actual response is 
correct and decreases for an incorrect actual response. As a consequence, a student likelihood 
will increase if all the actual responses are correct and decrease if all the actual responses are 
incorrect. This in turn pushes the location of the maximum likelihood solution for the given 
student to plus or minus infinity. For the item likelihood, a similar statement holds. For a generic 
response pattern, the student likelihood can be well-approximated by the density function of the 
normal distribution (especially when the number of items is large enough) and its curvature will 
be inversely proportional to the asymptotic standard error of the ability estimates. 


3 Multidimensional Item Response Theory 


3.1 Basic Models 

Even though widely investigated, MIRT is not yet widespread as an operational model. 
Hence, identifying the major players among the competing MIRT models is difficult. Here, only 
two models are discussed, one by Whitely (1980) and another one by Reckase (1997). With an 
item, a vector of discriminations a £ M* 0 and a vector of difficulties b £ are associated. With 
these the functional representation of the dimension-wise independent MIRT model of Whitely 
has the form 


D 


I a,b 


d D 


[o,i], 6»% b (0)= n 


i 


< 2=1 


1 _|_ e -a d (6 d -b d ) ' 


( 8 ) 
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If the conditional probability of passing the dth dimension of the item is given by 1+e _ ad ^ d _ bd) , 
then (8) can indeed be understood as the joint probability of passing all the independent 
dimensions of the item. Unless there are separate observed scores for each dimension, language 
such as correct response on dimension d cannot be used. To address this lack, this paper uses the 
phrase passing a dimension term , which may refer to an unobservable event. 

McKinley and Reckase (1982) (see also Reckase, 1997) put forward a model that takes the 
functional representation 


fSP . 
J a,b ’ 


r d 


[o,i], 


* - /«?*(*) 


i 

1 -|_ g— {fl | 0>“6’ 


(9) 


where a is as before and b G M. (x \ y) = Y^d=i x dUd is the usual scalar product of x, y E If 0 . This 
paper uses the term scalar product MIRT to refer to this model. See Figure 1. 



Figure 1 Dimension-wise independent and scalar product MIRT hypersurface. 


3.2 Definition of Multidimensional Item Response Theory 

The goal in this section is to give a definition of MIRT with as few assumption as possible. 
MIRT postulates that with a single item, multiple cognitive abilities could be detected. To 
accommodate this idea, one has to change the model for the ability space from the one-dimensional 
vector space R to a finite-dimensional vector space Vg. While any finite-dimensional vector space 
V is linearly isomorphic to for D = dirn(U) [see (10) for an explicit way of constructing such 
an isomorphism], this isomorphism is not canonical (there is no unique isomorphism V —> M D ). 
Because of this and other reasons that will become clear, is not used here as a mathematical 
model of ability space. 
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The reader unfamiliar with these notions is referred to Halrnos (1974) for an excellent 
introduction to linear algebra. Also, an intuitive understanding of the basic notions of smooth 
manifolds should help understanding of what follows, although it not strictly necessary. Among 
the many fine references to the topic, the interested reader may find Warner (1971) useful. 

The basic object in uni dimensional IRT is the item response function (IRF) and its 
graph, the item response curve (IRC). The graph of a function / : A —> B is a subset 
graph(/) = {(x, f(x)) £ A x B \ x £ A} of A x B. IRC is a one-dimensional smooth submanifold 
of i x [0,1]. While there is a scaling freedom even in the one-dimensional case (e.g., the 
[in]famous 1.7 multiplier in the logistic models), the possibility of ambiguous interpretation is 
minimal and one may use the functional (IRF) and the geometrical (IRC) representation almost 
interchangeably. 

In the multidimensional case, however, the matter is not so straightforward. As this paper 
shows, the functional representation and the geometric representations are different in a subtle 
way. One way to keep the presentation coordinate-free in multidimensional IRT is to postulate 
that the theory is given by an item response hypersurface (IRHS'). As in the unidimensional case, 
the IRHS is used to express the probability of correct response given an ability in Vq. 

Before this paper defined this notion, some notations should be fixed. For any v £ Vq, the ray 
of v is defined to be the line M • v in Vq determined by v: R • u = {Aw £ V<? | A £ R}. Similarly, for 
v,w £ Vq, the w-directed line going through w is defined by 

w + R • v = {u> + Xv £ Vq | A £ R}. 

For the notion of IRHS, there is then the following: 

Definition 1 A dichotomous item response hypersurface (IRHS) is a D = dim(VJj) dimensional 
smooth submanifold M of Vq x [0,1], so that for any two vectors v,w £ Vq the intersection of 
(w + R • v) x [0,1] and M is a graph of a monotonic function w + R • v —> [0,1]. 

A MIRT model is given when an IRHS is given. 

Note, that while w + R • v is not canonically isomorphic to R, nronotonicity of the map 

fv,w : w + R • v —► [0,1], Ah f V)W (w + Xv) 

can be unambiguously defined by requiring that either f VjW (w + Xv) < f v . ul (w + fiv) or 
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fv,w(w + Xv) > f v>w (w + gv) for all A, g E M whenever X < fi. The notation f v = f v> q is used here. 
Figure 2 shows the intersection of {w + M • v) x [0,1] and M in two dimensions. 

To understand the definition better, first assume that v is arbitrary and w = 0 in the 
definition above. Then, the line w + M • v = M • v can be understood as an ability direction. The 
monotonicity requirement of Definition 1 asks for the natural feature that, as the ability given by v 
increases, the probability of the correct response increases as well. For nonzero w, the requirement 
is equivalent to the conditional probability of correct response being monotonic with respect to 
one ability when the rest of the abilities are fixed to a certain, not necessarily zero value. To 
be precise, one should say that for w ^ 0 there exists a basis of Vg so that the monotonicity 
requirement reads as the interpretation above. Furthermore, for any basis of Vg, Definition 1 will 
ensure the monotonicity of the conditional probability of correct response for any ability direction 
given any fixed values for the rest of the ability directions (as defined by the basis). 



Figure 2 Intersection of (rc + M u) x [0,1] and IRHS. The intersection is the bold curve, 
which is required to be monotonic. 

Note also that the collection of maps f VjW for v,w S Vg defines the IRHS completely. For this 
reason, the notation / M , or / if no confusion may arise, is used for the function describing the 
IRHS M C V 0 x [0,1]. 
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One may be tempted to object to the use of notions like manifold and hypersurfaces. It is 
very important to note, however, that the conditional probability of correct response has been 
given by a hypersurface in the usual MIRT literature as well. One major difference in terminology 
is that it was still called surface in any dimension, which is a correct usage only in Dimension 2. 
In higher dimensions, the object at hand is a hypersurface, a special case of higher dimensional 
manifolds. 

A basis D = (iq,..., vd) in Vg defines a unique isomorphism 

D D 

iv '■ Vg ] A iVi e-> ] A iei, (A i G M), (10) 

i=l i= 1 

where (ei)f =1 is the standard basis of R d : (a)j = Sij (Sij is the Kronecker delta). This 
isomorphism can be trivially extended to a diffeomorphism 

■00 : Vo x [0,1] -*■ R D x [0,1], (v,t) ^ (i 0 (v),t) (11) 

and via this diffeomorphism the IRHS may be transferred from Vg x [0,1] to R D x [0,1]. Now, in 
R d x [0,1] the image of the IRHS may be given by the graph of a smooth function / : R D —> [0,1]. 
Note, however, the important difference between using a functional representation like this latter 
one and using the hypersurface representation directly in Vg x [0,1]. The functional representation 
depends on the basis chosen to establish the diffeomorphism ip 0 and different bases may result in 
different functional representations. 

It is tempting to extend this definition to polytomous multidimensional items by defining the 
polytomous collection of item response hypersurfaces for a polytomous item by requiring that the 
above discussed intersection be a collection of unidimensional polytomous item response curves 
as produced by some unidimensional polytomous IRT model (e.g., Muraki’s [1992) partial credit 
model). The investigation of this possibility is postponed for a forthcoming paper. 

3.3 Properties of Item Response Hypersurfaces 

This section proves the main theorem of the paper. For the sake of transparency, this section 
starts with the two-dimensional case, which is then followed by the more involved general theory. 

3.3.1 Two-Dimensional Case 

Using the nronotonicity of the model, it is possible to prove an interesting elementary 
property. 
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Lemma 1 In any two-dimensional MIRT model, there exists a line in Vg through the origin so 
that f v is constant. 


Proof: Choose a vector, v £ Vg. Note that if inf f v (Xv ) = sup f v (Xv), the lemma is proved and 

asm 

the sought after line is M • v. Therefore, one may assume that inf f v (Xv ) < sup f v (Xv). For such a 

A SR AeR 

vector, either lim f v (Xv) = sup f v (Xv) or lim f v (Xv) = sup f v (An). Let P be the set of vectors 
A^oo AeR A^-oo AeR 

satisfying the first and N be the set of vectors satisfying the second condition. Both of these sets 
are nonempty, and by continuity, both sets are open. Also, they are clearly disjointed. Therefore, 
there is a vector u £ Vg so that u ^ N U P. Along the line M • u the function / is constant. ■ 

Note that the proof only uses nronotonicity with w = 0. Utilizing it for general w, the same 
argument provides the following: 

Lemma 2 In any two-dimensional MIRT model, through any point w £ Vg there exists v £ Vg so 
that, along the v-directed line going through w, the function f v , w is constant. 

The term w-constant line, or simply constant line, is introduced here for the ^-directed line going 
through w as in Lemma 2. 


f=const. 



Figure 3 Nonunique constant line results in constant MIRT model. 

Analyzing the properties of these constant lines further, one can see that they are actually 


parallel to one another, which leads to the following: 



Lemma 3 Let w, w' £ Vq be two points. Let v, v' £ Vq be the corresponding directions of the two 
constant lines. Then v = pv' for some p £ M. 

Proof: First, note that if there is a point w £ Vq so that there exist two rc-constant lines, then the 
model is trivial (/ is constant) and the statement is true. For, let w' and w" be the intersections 
of a general position line in Vq with the two re-constant lines, respectively. Because / is monotonic 
along this line, f(w') = f(w"), f is constant between w' and w", that is, f[tw' + (\ — t)w") = f{w') 
for all t £ [0,1]. Using this argument for every line in a general position proves that / is constant 
everywhere (see Figure 3). 

Now assume that the constant lines through w and w' are unique. If the two lines are not 
parallel, then they will have an intersection, and an argument similar to the previous one shows 
that / is constant. ■ 

The corollary of the previous observation is given in Theorem 1. 

Theorem 1 Any two-dimensional MIRT model is a trivial extension of a unidimensional IRT 
model. 

Proof: Lemma 3 showed that a two-dimensional IRHS is nothing but a collection of parallel lines. 
Let v £ Vq be the direction of these lines. By choosing a transversal M • u (a line that intersects 
all of them) to this collection, the IRHS can be given by the function f u : M • u —> [0,1]. One can 
express an arbitrary w £ Vq as a unique linear combination w = pu + Xv and write 

f(w) = f{pu + Xv) = fu^pu)- (12) 

This function f u can be thought of as a unidimensional IRT model. ■ 

3.3.2 D-Dimensional Case 

Technically, the D-dimensional case is not that much more complicated than the two- 
dimensional one. It is just much more difficult to visualize the corresponding geometric objects. 
As pointed out earlier, the conditional probability surface is not two-dimensional, so strictly 
speaking it is not a surface in higher dimensions. Three-dimensional training does not allow one 
to “see” objects in higher dimensions. The formalism built in the previous section, however, is 
applicable, with appropriate modifications, to this situation as well. 
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The proof of Lemma 1 works for any dimensions. Applying the nronotonicity argument for 
arbitrary (v, w) as above proves the corresponding 

Lemma 4 In any MIRT model there exists a hyperplane H w in Vg through w £ Vg so that fu w is 
constant. 

Proof: Here, fu w is the restriction of / to the hyperplane H w . As before, w = 0 is proven explicitly; 
the general case follows the same argument. As before, the open sets P and N are denoted, with 
P = —N. Exclude the trivial case of P = 0. It is clear that PUlV ^ Vg. Locally the boundary of 
P (the closure of P minus P) is a D — 1 dimensional submanifold (D = dirnVe). Therefore there 
exists a collection (ci,..., cd~\) of points in Vg\(P U N) so that (ci — w ,..., cd -i — w) spans a 
hyperplane H w . Now, along the line segment joining any point on M • (cj — w) with any point on 
M • (Cj — w) for some i / j, the restriction of / should be constant (Figure 3). Repeating this 
argument for each pair of line segments shows that along the entire hyperplane / is constant. If 
there is another hyperplane with this property, then P = 0, which is excluded. ■ 

Now, a D — 1 dimensional hyperplane is to a D-dimensional space as a line is to the plane. 
Using this intuition, it is not difficult to adapt the formal proof of Lemma 3 to prove Lemma 5: 

Lemma 5 Let w,w' £ Vg be two points. Let H w ,H w i C Vg be the corresponding two constant 
hyperplanes. Then H w and H w i are parallel. 

Now one can rephrase the main theorem in arbitrary dimension. 

Theorem 2 Any MIRT model is a trivial extension of a unidimensional IRT model. 

For the sake of explicitness, write f M for an arbitrary IRHS in terms of univariate IRT model. A 
transversal u £ Vg is fixed to the collection of constant hyperplanes. First, observe that for any 
w £ Vg there is a unique decomposition w = /m + Xv with v £ H w . Then 

f M (nu+\v) = fu(iiu). (13) 

Note that if the usual 2PL or 3PL models are chosen, the construction yields the scalar 
product model. It is also interesting to note that the MIRT generalization of the Rasch model is 
equivalent to the generalization of the 2PL model. This is because of the following: while within 
the univariate Rasch model, one may assume that the slope is fixed, but when more dimensions 
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are considered simultaneously, the assumption of equal slopes is not valid. The relative positions 
of slopes to one another should be determined during the estimation procedure in lack of a priori 
information. 

This kind of models were called generalized compensatory models (GMIRT) in Zhang and 
Stout (1999). The link function of an IRHS as GMIRT is fff . 

3.3.3 Absolute Functional Representation for the Scalar Product Model 

A notable feature of the scalar product model is that using the dual of a vector space it can 
be defined without referring to coordinates even in its functional form. First, recall that the 
dual V* of a finite-dimensional vector space V is the finite-dimensional vector space of the same 
dimension of linear maps V —> M : 


R* := {p : R —> M | p is linear}. 


(14) 


The duality is the obvious map 


( | ) : V* X V —> M, (p, v) i—> ( p | u) := p(v). 


(15) 


That is, for any p £ V* and v G V the quantity ( p \ v) is a real number. It is important to note 
that the duality, unlike a scalar product, does not involve any choice. 

Now, if in MIRT the choice is made that the ability is modeled by the vector space Vg as 
before and the item is modeled by the discrimination a € Vf in the dual space and a real number 
b, then the IRHS of the model is given as the graph of the following function: 


fU ■ Ve 


[o,i], fim : = 


1 


1 -I- e -(a I 9)-b' 


(16) 


In addition to its very satisfying and elegant nature, this model has the computational advantage 
of having the same functional representation in any coordinate system. As is shown later, the 
dimension-wise independent model does not share this nice invariance property. 


3.3.4 Interpretation of the Main Theorem 

The statement of the main theorem excludes many existing MIRT models from the pool of 
monotonic MIRT models. The author’s reading of the main theorem is that the only relevant 
MIRT model is the one defined in (13). This interpretation is backed by the fact the useful 
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estimation methods exist only for the scalar product model, the most relevant of the above 
extensions (Reckase, 1997). In the view of Theorem 2, there seems to be a good reason behind 
that. It seems that lack of monotonicity prevents one from maximizing the likelihood function of 
MIRT models excluded by this approach. This certainly defines a valid future research direction. 
Also, the existence of an elegant coordinate-free functional representation makes the scalar 
product model even more appealing. 

On the other hand, model building always has many steps that cannot be entirely backed by 
theoretical considerations. The process sometimes is dictated by personal preferences and tastes. 
It is possible that some readers may not be willing to accept the requirement of monotonicity 
as formulated in Definition 1 as a crucial and necessary feature of an MIRT model. For those 
readers the main theorem is interpreted a bit differently. First, note the close connection between 
the notion of compensatory model to monotonicity. Usual terminology is that the model is 
compensatory if the probability of the correct response may be high even with the lack of ability in 
all but one dimension. That is, sufficiently high ability in one dimension is able to compensate for 
the lack of it in other dimensions. In fact, compensatory property follows from monotonicity as an 
easy application of Theorem 2. If compensatory property is understood in a sense that it is true 
in any coordinate system, then the reverse is also true, and the two notions are equivalent. With 
this in mind the theorem states that any compensatory MIRT model is a direct generalization of a 
univariate IRT model. 

In either way, Theorem 2 establishes a prominent role for the scalar product model as an 
MIRT model. 

3.4 Estimation in Multidimensional Item Response Theory 

The paper now restricts attention to the scalar product model. A typical student likelihood 
(with D = 2) is given as the right-hand side graph in Figure 4. As in unidimensional IRT, the 
maximum place of this function plays a special role in the estimation of MIRT model parameters. 
A curious feature of this graph is that a pronounced imbalance can be observed between the 
standard errors of the two ability estimates. Here, standard error is understood as the inverse of 
the curvature of the graph at the maximum place. There is a well-identified direction in which 
the standard error is minimal, and in the direction orthogonal to this, the standard error appears 
to be much bigger. One may even say that, despite these efforts, the model shows definite signs 
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of unidimensionality. The reason behind this is very simple. A student likelihood is formed as a 
product of probabilities of the actual responses given by item response hypersurfaces similar to 
the one shown on the RHS of Figure 1. These hypersurfaces are always increasing towards the 
first quadrant (correct response) or towards the third quadrant (incorrect response). Hence, the 
product of these will be the ridge in Figure 4. It is a ridge because the observed response is either 
correct or incorrect, and no distinction is made between events of the students using only one of 
the dimensions correctly during the assessment. In other words, since there is no observed data for 
the different dimensions, the model will not be able to provide two distinct, meaningful estimates 
for the abilities of the person on the different cognitive dimensions. 



Figure 4 Scalar product MIRT student likelihood. 

3.5 Dimension- Wise Independence 

The careful reader should have noticed that, concerning one particular point, the presentation 
is not faithful to its own principles. That is, the notion of dimension-wise independence was used 
without any discussion of its invariance or coordinate-system independence. It is easy to see that 
the dimension-wise independent model does not satisfy the requirement of nronotonicity, therefore 
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it would not be considered it as a valid MIRT model. On the other hand, it might be useful to see 
explicitly how badly the the functional representation of the dimension-wise independent model 
behaves and so appreciate the niceties of the scalar product model even more. 


Invariance of dimension-wise independence for the the model 

k„ ■■ r b - [o, i], ftm = n 


i 


d= 1 


l _|_ e -a d (6 d -b d ) 


(17) 


would mean that the factorization property holds in any other coordinate systems. 

Mathematically, this would require that, for any invertible matrix G £ GL(D) (G expresses 
change of coordinates), is a function 

h G : M x M x M —M, (<xi, 6i, t) i—>• h G {a\,b\,t) 

with a pair of invertible matrices U, V £ GL(D), so that when 9 = G ■ 9' (9, 9' £ M 15 ), there is have 
a factorization 

D 

f:, b (9) = Hh G «Ua) d ,(Vb) d ,9' d ), (18) 


d= 1 


that is, 


= faAG- 6 ') 

D 

- n 


L = ^ 1 -g e ~ a d(T, d ' = 1 9dd' e ' d t- b d) 
D 


Ylh G (a' d ,b' d ,9 ' d ), 


(19) 


d= 1 


with a' d = ( Ua) d and b' d = (Vb) d . The role of U and V is to ensure that the function h G is 
the same for all factors in the product by allowing this function to depend on different linear 
combinations of the elements of a and of b. 

To show that this is too much to ask for in general, first assume that a factorization 
f(x, y) = h(x)g(y) holds for some function / so that h( 0) / 0 and g( 0) / 0. Then, 

/(*, 0 ) 


and 


h(x) = 

g(y) = 

1 


9(0) ’ 
f(0,y) 

HO) ’ 

_ = f(x,y) 

h(0)g( 0) f(x, 0)/(0, y) 


( 20 ) 
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Now, for the sake of concreteness, take D = 2 and a = (ai,ai) £ M 2 and b = (0,0) £ M 2 . Also, 
1 1 


take G = 


1 -1 


. With these, (19) becomes 


1 


with some h, g 


1 + e a i( 0, i+ e 2 ) 1 + e 

From (20) the function 

1 


-a l( e'-e') = )M) 


( 21 ) 


l_l_ e a i( 0 i+ 0 2^ 1 +e ai ^i e 2 ) 


/i(0)fl(0) 


( 22 ) 


(l+e a i 0 i ) 2 1 +e “ lS 2 i+ e + a l®2 


should be constant. This is clearly not the case, showing that the factorization (19) does not hold 
in general. 

It seems that the definition of dimension-wise independence is not an absolute one. One can 
either drop it altogether, or if need arises, change it to the following: 


Definition 2 An MIRT model given by an IRHS is dimension-wise independent if there exists a 
coordinatization of abilities so that the functional representation of the model f a ,b(6) can be written 
as a product of factors 

D 

fa,b(8) = l[h(a d ,b d ,e' d ). (23) 

< 2=1 

The specialty of this property comes from the fact that for a general IRHS it is very rare 
that the functional representation can be factored so that one may consider it dimension-wise 
independent. This interpretation was used throughout the paper, when the Whitley model was 
called dimension-wise independent. 


4 Conclusion 

A coordinate-free definition of MIRT has been put forward in the paper. The main argument 
is that in a coordinate-free setup it is easier to distinguish genuine MIRT objects from potential 
artifacts. These artifacts can be notions and relationships, either of which should not be considered 
integral parts of the model since key features apparent in one coordinate system could vanish in 
another. This paper showed that it is possible to provide a full classification of monotonic models 
solely based on general, coordinate-free considerations. 

It is very important that the reader does not mistake the promotion of the coordinate-free 
description as an argument for a completely coordinate-free handling of the entirety of MIRT. In 


15 



fact, meaningful MIRT practice cannot exist without a choice of coordinates. In addition to this, 
every discussion of MIRT features can be fully carried out using W D as the main model space for 
abilities. Should such a path be chosen, however, one has to be careful to meticulously maintain 
the coordinate-system invariance of the theory every step of the way. The contribution of this 
paper is an introduction of a framework to ease this burden by keeping the presentation absolute 
(without choosing any coordinates) for as long as possible. The paper shows that one may be 
able to formulate general statements and reach valuable insights before switching to relative mode 
by an introduction of a particular basis. It is likely that someone may observe the relevance of a 
notion while using a particular coordinate system and may want to establish whether the notion 
is invariant by trying to create a definition in the absolute framework presented here. 

It is noteworthy that the necessity of the existence of a coordinate-free representation of our 
physical world led Einstein to formulate both the special and the general theories of relativity 
(1905, 1916). The fundamental dogma in relativity theory is that the events of the physical world 
take place without there being aware of any coordinate system. Therefore, any faithful description 
should be invariant of the change of coordinate system. Better yet, a description of the physical 
world is sought that bypasses the use of coordinates altogether. 

A reader interested in the successes of coordinate-free description of the physical world may 
also find the books by Matalocsi (1986, 1993) useful. 
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